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Tropospheric inhomogeneities limit the accuracy with which a path delay in an 
arbitrary direction can be estimated from calibration measurements in different 
directions. This article demonstrates a mathematica/ procedure that has the po- 
tential for minimizing errors in the estimated geometrical and tropospheric path 
delays. The error is minimized by applying least-squares estimation to a combined 
set of observables in the calibration directions and the direction to be calibrated. 
A simulated test of this procedure was conducted using a model set of error- free 
calibration measurements. In the absence of geometrical delay mismodeling, the 
simulation yielded delay errors which vary from about 1 mm at zenith to about 
1 cm at 10 deg. The main principles of how this procedure could be applied to 
improve accuracy of deep space tracking using Global Positioning System (GPS) 
data are also discussed. 


I. Introduction 

Uncertainties in tropospheric path delays are a major 
source of error in deep space tracking. Inhomogeneities in 
tropospheric water vapor can result in zenith path delay 
calibration errors at about the 1-cm level, and inhomo- 
geneities in the dry troposphere at about the 1- to 3-mm 
level, over a period of several hours [1]. The 1-cm error 
limits the tracking accuracy of DSN-based very long base- 
line interferometry (VLBI) to about 1 nrad for the an- 
gular position (at the intercontinental baselines of about 
10,000 km) and to 2 x 10“ sec/sec for the delay rate 
(at zenith) for a 1000-sec scan [1]. Future missions would 
benefit from troposphere calibration at the 1-mm level. 

The error in the estimated path delay is determined 
by a variety of error sources whose relative importance 
depends on the calibration instrument. For example, for 


instruments which measure radio emission, such as water 
vapor radiometers (WVRs), the error is limited by the 
accuracy with which path integrals involving the imagi- 
nary part of the index of refraction can be related to in- 
tegrals which involve the real part of the refraction index 
[2]. Whereas WVRs may be pointed in the direction to 
be calibrated, this may not be true in general for other 
instruments. If off-1 ine-of-sight measurements are used in 
the calibration process, the error will inevitably be affected 
by tropospheric inhomogeneities, and its magnitude will 
depend on mathematical analysis of calibration data. 

This article suggests and illustrates a mathematical 
procedure that minimizes the tropospheric inhomogeneity- 
induced error for path delays inferred from observables in 
other directions. The procedure involves the application of 
least-squares estimation to a combined set of observables 
in the calibration directions and the direction to be cali- 
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Fig. 5. The phase time series of Fig. 3 with HeNe LDI phases, scaled to infrared cycles, 
differenced between telescopes, and multiplied by 3. 
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Fig. 6. The calculated correlation between telescope-differenced 
HeNe LDI data and interferometric data for the conditions of the 
data in Figs. 3-5. The horizontal line Is the actual correlation from 
the data of Fig. 5. 
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Fig. 7. The calculated astrometric accuracy, as a function of scan 
length, with and without optimal HeNe LDI calibration, for the con- 
ditions of Figs. 3-5. 










brated, and the use of the observable variance-covariance 
matrix during estimation. In order to demonstrate the 
procedure simply, this article makes the following assump- 
tions: The calibration instrument is error-free (i.e., no 
other error source besides tropospheric inhomogeneities is 
present), and it detects path delays. These assumptions 
are used for clarity of the presentation; the method may 
be generalized to other types of measurements. An exam- 
ple of how it can be adopted to more realistic (noisy) data 
is discussed in Section IV of this article. If all other er- 
ror sources are neglected, the calculated error determines 
the maximum achievable accuracy for path delay estimates 
using observables in other directions. Another (more com- 
monly used) procedure determines the delay in the zenith 
direction by averaging over delays projected into the zenith 
from many calibration directions. This technique will be 
referred to in this article as “zenith mapping,” Strictly 
speaking, zenith mapping produces minimal error only 
for a horizontally homogeneous troposphere. This article 
will quantify the error reduction obtained from the use of 
the present technique for an inhomogeneous troposphere. 
Since the technique’s principal intended application is to 
improve the accuracy of deep space tracking, the direction 
to be calibrated will be sometimes referred to in this article 
as the DSN antenna pointing direction. 

Section II describes the mathematical procedure used 
to minimize the estimated delay error. The procedure is 
exemplified in Section III for a model set of calibration 
directions (assuming error-free instrumentation) that co- 
incides with directions of lines between a ground-based re- 
ceiver and Global Positioning System (GPS) satellites visi- 
ble at Goldstone.^ These directions were selected as a mat- 
ter of convenience and because they are known with great 
accuracy.^ Tropospheric inhomogeneities were assumed to 
be generated by Kolmogorov turbulence and transported 
past the observer by the wind [1]. Section IV discusses the 
possibility of adapting the procedure to GPS data. Sec- 
tion V is a summary with recommendations for further 
studies. 

II. Mathematical Approach 

This section outlines the main principles of the mathe- 
matical procedure that minimizes the estimated delay er- 


^ Listings of GPS coordinates were provided by G. Purcell, Track- 
ing and System Application Section, Jet Propulsion Laboratory, 
Pasadena, California, October 1992. 

^ Other off-line-of-sight measurements to which the method may be 
^ldapted include microwave temperature profiler [3], lidar [4], radio 
£icoustic sounding [5], and tracking using a number of proposed 
communication satellites systems, such as Motorola’s IRIDIUM 
constellation. 


ror in an arbitrary direction. The main principle of the 
procedure is the application of least-squares estimation to 
a combined set of observables in the calibration directions 
and the direction to be calibrated (also referred to as the 
pointing direction of the DSN antenna), and the use of the 
covariance- variance matrix of observables to weight the 
quadratic form of observable residuals during estimation. 
To simplify the illustration, all measurements are assumed 
to be error-free and to have produced path delays. 

The coordinate system is shown in Fig. 1. A calibration 
instrument, located at the distance R from the axis of the 
DSN antenna, measures tropospheric path delays in 
N different directions (Ei, (f>i)j where Ei and <^i are eleva- 
tions and azimuths, respectively, and i = 1,..., V. A DSN 
antenna (the direction to be calibrated) points in the di- 
rection of the elevation, E*,, and azimuth, <f>g^ By assum- 
ing that the ionospheric delay has been calibrated (e.g., 
by using two frequency measurements), the DSN antenna- 
measured delay is Ltot,s = Lg^s + where 

and Lir s are the total, geometric, and tropospheric delays, 
respectively. The combined set of observables in the DSN 
antenna and calibration directions is; 

^tot,s — ^tr,s “H (^^) 

Li = Ltr^i i=I,...,V (lb) 

where the symbol Li designates the observable Ltj- i. The 
delay of interest for deep space tracking is the geometrical 
delay j. To obtain the best estimate Lg^^^ we (1) note 
that Lir^s and Ltr/s are related through tropospheric cor- 
relations, (2) parameterize Eqs. (la) and (lb) with the 
help of the statistically averaged zenith delay Ltr^z (that 
is, is not the instantaneous zenith delay, but rather 
a delay averaged over all possible tropospheric patterns 
for the site [1]), and (3) apply least-squares estimation to 
the parameterized Eqs. (la) and (lb). Assuming that 
long-range refraction gradients [6] are absent,^ the param- 
eterized Ltr,s — T ^tr^s ^-Ud Li^ j — i^Lif i) -|- 

where {Ltr,s) = A^Ltr^z and {Ltr,i) = AiLtr^z, ( ••• ) des- 
ignates the expectation value, Ag and Ai are air masses, 
and and represent tropospheric inhomogeneities 
in [Eg^<j)g) and [Ei,4>i) directions, respectively. Note that 
because of (2) above, correlations of and ^tr,i can be 
evaluated in the statistical sense. By defining the observ- 
able, parameter, and tropospheric inhomogeneity column 
vectors F = [Ltot,JA,, Li/Ai\, X = [Ltr.z, 


^ Unmodeled horizontal gradients may be on the 1-cm level for dry 
delays at 10-deg elevation [6]; the error can be reduced by gradient 
modeling during estimation. 
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and e = [cfr^s, respectively, and by assuming that 

€ir/^ have zero means and a variance-covariance matrix 
the parametrized Eqs. (la) and (lb) are solved by 
minimizing the quadratic form of the weighted residuals 
{F — AXY W {F — AX), where A is the mapping matrix 
m F ~ A X e, 


1 0 

\1 0 / 


where c; = ~{A,/Ai) {{A'^WA)~^A'^W)2,i+\ represents 
contributions from the zth calibration direction. By us- 
ing Eq. (5), the estimated delay error can be expressed as 
follows: 


2 _ 2 
L,tr,s L,g,s 




N 


The result is [7] 

X = {A^WA)-^A^W F (3) 


N 

^tr,j ) ( 6 ) 


where X is the column vector X = [Ltr,zy and 

the superscript T designates the transpose matrix. The 
error in is given by the square root of the matrix 

element (< 7 ^) 2, 2 of the variance cr^ \ 

al = ({X-X){X-Xf) 

= B-^A'^W COW {F,F'^) W A (4) 


where cov (F, F^) is the observable covariance- variance 
matrix, and B = A'^ WA^ 


The above-described procedure yields the best estimate 
and Ltr^z- Note, however, that once Lg^g has been es- 
timated, Eq. (la) can also be used to estimate the actual 
line-of-sight tropospheric delay Ltr^g. By using Eq. (la), 
the estimate Ltr^s ^ ^tot,s — Ltr.s is the best 

estimate for Lir^s because Lg^g is the best estimate for 
Lg^g, and because the measured Liot^s contains the effect of 
line-of-sight inhomogeneities (note that Ltr^s differs from 
Lir^zAg, which is the best estimate for the statistically av- 
eraged delay). Note also that in the absence of other error 
sources, the error in Ltr^g is equal to the error in Lg g. 


It will be useful for the ensuing discussion to explicitly 
write Lir^s as"^ 


N 

Ltr^s = ^tot^s ~ ^toi^s AgX2 — ^ ^ (b) 

1 = 1 


Equation (6) displays an explicit dependence on the 
coefficients c^’s and on correlations between tropospheric 
path delays. When c^^s are determined by using some as- 
sumed W~^ (the so-called consider analysis [1]), the error 
will be bigger than the error obtained by using the ac- 
tual observable [7]. A W~^ often used during the 

estimation procedure is the diagonal matrix, W~j^ = 6ij 
(where is the Kronecker delta). The use of a diagonal 
W~^ corresponds to assuming that observable errors are 
uncorrelated, and Eq. (3) yields^ c, = Ag/{NAi), which 
are the same Cj’s as those which would be obtained if one 
set ^tr,z = ^<r, i/(Ai N), and mapped Ltr^z from the 

zenith to the {Eg,<j>g) direction by using the air mass Ag. 
The use of zenith mapping will minimize the error only 
for a horizontally homogeneous troposphere. For an inho- 
mogeneous troposphere, the error is minimized by using 
the observable W~^ = cov(F, F^); c^s will then depend 
on the full (in elevation and azimuth) angular separations 
between the observed and calibration directions (includ- 
ing the offset R). By setting W~^ = cov (F, F^), Eq. (4) 
reduces to = {A^W A) ^ 


III. Results for an Error-Free Calibration 
Instrument 

Equation (6) was evaluated for a model set of calibra- 
tion directions assumed to coincide with directions of lines 
between a ground-based receiver and GPS satellites visible 
at Goldstone. These directions were selected because there 
are between 6 and 10 satellites visible from any ground- 
based site, and the satellite trajectories are known with 


^ The proportionality coefficient between Lg^g and Ltot,s is equal to 
1 because no other observable depends on Lg^a> 


^ For elevation-independent correlations, W- ^ = 6tj/AtAj, the co- 
efficients would become c, = AgAij 
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great accuracy. Other directions could be selected; con- 
clusions similar to those derived here will apply to all 
other selections. The error was quantified by using the 
Kolmogorov turbulence model for the evaluation of cor- 
relations between wet troposphere inhomogeneities, with 
numerical constants given in [1] (see also Appendix A of 
this article). The inhomogeneities were frozen into the 
troposphere slab transported p 2 ist the observer with the 
wind velocity u = 10 m/sec. The structure function con- 
stant C — 1.1 X 10“ corresponds to average DSN 
observation conditions of about a 6-cm zenith wet delay 
for = 2-km- thick wet troposphere slab [8]; in more hu- 
mid weather, the structure constant will be greater, de- 
pending on turbulence. The dry fluctuation was assumed 
to be 30 percent of the wet fluctuation [1], which for un- 
correlated fluctuations contributes less than 10 percent of 
the total error for correlations added in quadrature. The 
error evaluated with the optimized Cj’s was compared to 
the zenith mapping error; the zenith mapping error was 
greater, especially at low elevations, as expected. The er- 
ror was quantified for a number of directions <^,), for 
satellite constellations stepped by 6 min during an 8-hr 
period following 12:00 a.m., July 23, 1992, and for four 
values of the separation = 0, 200, 500, and 1000 m. The 
biggest contribution to the estimated delay comes from a 
calibration direction nearest to the direction to be cali- 
brated; the error vanishes when the two lines coincide. In 
the absence of a clearly nearest direction, all calibration 
lines contribute, increasing the error. 

A. Delay Error 

The error was found to depend in a relatively well- 
defined manner on the elevation , but not (because the 
satellites are distributed over the azimuth fairly uniformly) 
on (f>s nor on satellite constellation. Figures 2(a) and 2(b) 
plot the azimuthally averaged error (dashed line) versus 
Eg for = 0 and 1000 m. The error spread due to (1) a 
360-deg range in 0,, (2) a 360-deg range in the azimuth 
of R, and (3) changes in the satellite constellation is also 
shown. The spread’s lower bound is 0 for all Eg in Fig. 2(a) 
{Eg extends from 10 deg in Fig. 2(a) because the satel- 
lites were cut off at 10 deg). The bound is 0 because at 
R = Qj any direction Eg will, for some combination of 
<j)g and satellite constellation, eventually coincide with one 
of the calibration lines, in which event the tropospheric 
inhomogeneity-induced error vanishes. Note that the av- 
erage error is closer to the error upper bound than to the 
lower bound, indicating that the error is nearly equal to 
the error upper bound in the majority of observations. 
Comparing Figs. 2(a) and 2(b), the average error has in- 
creased only very little, while the lower bound has become 
nonzero by increasing R from 0 to 1 km. The average 


error has increased little because the DSN-observed and 
calibration lines intersect (and thus probe the same sky 
region) even when = 1 km; the lower bound has become 
nonzero because there are no coincident lines when i? > 0. 
The average error is fairly flat (about 1 mm) for all Eg 
between zenith and 40 deg, and it increases rapidly with 
decreasing Eg at lower elevations; the error is about 1.6 
and 12 mm at Eg ~ 30 and 10 deg, respectively. 

Figure 3 illustrates how close the calibration and DSN 
line must be for the error to be less than some desired 
value. The error (the solid-line curve) will be less than 
0.5 mm when the angular separation between the DSN 
and the nearest calibration line (jR — 0 in Fig. 3) is less 
than 3 deg. How small the angular separation must be de- 
pends on elevation: modeling results suggest that for the 
error to be less than 1 mm, the separation must be less 
than 5 and 2 deg when Eg — 30 and 20 deg, respectively. 
Many estimation strategies estimate zenith delays by aver- 
aging over all calibration directions. The zenith mapping 
error (the dash-dot line) is bigger than the optimized, az- 
imuthally averaged error (the dashed line) by an amount 
which decreases monotonically with increasing elevation, 
until, near zenith, the errors are nearly the same (because 
the mapping distance is short there). 

Signal integration averages out the tropospheric in- 
homogeneity-induced error as the inhomogeneities are car- 
ried by wind. Figure 4 shows the effect of the signal in- 
tegration time T and wind velocity v on the error.® The 
error decreases and its spread (due to different wind direc- 
tions) increases slowly with increasing T. Note that the 
average integration time, Ti/ 2 j required to reduce the error 
to one-half of its instantaneous value is (for v = 10 m/sec) 
less than 8 min for all Eg > 30 deg, and it increases to 
about 12 min at Eg — 10 deg. To minimize the error, c,’s 
used for Fig. 4 were optimized for the wind. A similar cal- 
culation using Cj’s optimized for zero wind has produced 
curves (not shown) that look the same as those in Fig. 4 
except that the error was approximately 10 to 30 percent 
higher (depending on elevation). This relatively small in- 
crease is good news, since the determination and inclusion 
of wind distribution in the estimation procedure could be 
nontrivial. 

The use of more than one calibration instrument will 
(in principle) increase the probability that one calibration 
line will be close to the observed line. Figure 5 shows the 
minimum error for three instruments positioned in corners 


® H = 200 m is used in most figures in this article, since the az- 
imuthally averaged errors do not differ too much for all R*s < 
1 km, and since mounting a calibration instrument in the center of 
the DSN antenna is nontrivial. 
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of an equilateral triangle around the DSN antenna. Com- 
pared to one instrument (positioned at /Z = 0 and 200 m), 
the error is smaller by up to 50 percent at high elevations 
(where it is already small), but only by several percent at 
low elevations (where it would be needed the most because 
the error is big there). The result suggests that multiple 
on-site instruments will be of limited use. 

B. Delay Rate Error 

Figures 6 and 7 show the inhomogeneity-induced delay 
rate error. Rate measurements are used for navigation 
and gravitational wave searches; it is desirable that the 
inhomogeneity-induced rate error not dominate the total 
error on any time scale. The rate error was evaluated as 
the square root of the variance cr|(T 5 c), 

a|(T,e) = ((s(T,e)-S(r,e))'^ (7) 

where B(T,c) and B{Tsc) are the estimated and actual de- 
lay rates, respectively, over the scan duration Tgc (see [1] 
and Appendix B). The variance was evaluated by using a 
linear fit to three equally spaced points^ '^sc^ Fig- 

ure 6 shows the error evaluated by using the optimized 
c,’s, and compares it with the uncalibrated error and the 
error evaluated by using zenith mapping for two elevations, 
Eg = 10 and 60 deg. All three errors decrease with increas- 
ing Tgc (as they should). However, whereas the optimized 
error is actually slightly bigger than the uncalibrated and 
zenith mapping errors for extremely short scans {Tgc < 
10 sec), the optimized error decreases with increasing Tgc 
more rapidly than the other errors, becoming smaller at 
a Tgc that depends on Eg. For example, at Eg — 10 and 
60 deg, the error becomes smaller at Tgc = 500 and 30 
sec, respectively. For the error to be reduced by using the 
optimized q’s, the scan must be longer than some criti- 
cal Tgc, which is longer at low Eg than at high Eg. At 
Tgc — 1000 sec {R = 200 in), the average error is about 
4 X 10“^^ sec/sec and 7 x 10“^"^ sec/sec at Eg = 60 and 10 
deg, respectively. 

The rate error can also become very small when the 
angular separation between the observed and nearest cali- 
bration line is sufficiently small; however, the condition for 
the separation smallness is tighter than that for the delay 
error. This is illustrated in Fig. 7 for R — 0. For the error 
to be less than 10“^^ sec/sec, the angular separation must 
be less than about 0.5 deg when Tg — 1000 sec and less 
than about 0.05 deg when Tg = 100 sec. 

^ It can be shown that as the number of fitted points increases to 
infinity, the rate error changes by only about 20 percent [1]. 


IV. Application to GPS 

This section discusses how the mathematical procedure 
described in Section III could be adapted to path delay es- 
timation using GPS data. GPS tracking depends directly 
on path delays.^ However, because of the complexities of 
GPS delay modeling and data reduction, the present ap- 
proach to GPS data analyses, driven by the requirement 
to eliminate or estimate clock errors and geometrical pa- 
rameter uncertainties,^ entails estimation of instantaneous 
zenith delays [9]. Parameter estimation is performed with- 
out the use of the tropospheric covariance- variance matrix. 
The zenith tropospheric delay is constructed by spatial av- 
eraging over all satellites. This zenith delay could be used 
to estimate tropospheric delays in the direction of the DSN 
antenna by scaling it (using the air mass) into the DSN 
antenna direction. This article suggests that in the ab- 
sence of other error sources, line-of-sight delays could be 
estimated more accurately by modifying the estimation 
strategy to include the tropospheric correlations between 
different lines of sight and by applying the estimation pro- 
cedure to the combined GPS and DSN data. 

The coordinate system is the same as in Fig. 1 ex- 
cept that receivers (which are a part of a large global 
network consisting of about 40 continuously operating re- 
ceivers set up to estimate the geometrical parameters and 
to determine the clock uncertainties with the best possi- 
ble accuracy) are separated by large (several hundred-km) 
distances. For M receivers, one can introduce M sta- 
tistically averaged, mutually uncorrelated zenith delays, 
Ltrj^zi where j = The DSN antenna is assumed 

to be at the site i = 1. At its completion, the GPS con- 
stellation will include 24 navigation satellites at about a 
20,000-km altitude and equally spaced in six orbit planes. 
A receiver at any ground-based site will track between 6 
and 10 satellites distributed more or less uniformly in az- 
imuth and typically above 10 deg in elevation. The col- 
lected data set covers a time period long enough for the 
number of data in the set to exceed the number of solve- 
for parameters. Including the delay for the radio signal 
received by the DSN antenna, the observables are 

LtotfS ^ ^tr,5 F (^^) 


® GPS receivers are used in an automated operating mode, work 
in all weather conditions, and sense the total (wet and dry) path 
delays. The satellites transmit carrier signals at two L-band fre- 
quencies (1.227 GHZ €md 1 .575 GHz), so that the ionospheric delay 
can be calibrated. 

^ S. M. Lichten, personal communication, Tracking Systems and Ap- 
plications Section, Jet Propulsion Laboratory, Pasadena, Califor- 
nia, October 1992. 
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i =1,..., AT j=l,.., M 


(8b) 


which is sufficiently small to argue that the use of Cj’s opti- 
mized for zero additional uncertainty will not significantly 
increase the total error. 


where N is the total number (assumed to be the same for 
all receivers) of calibration data at one receiver site (the 
subindex i counts both the visible satellites and the data 
sequence for each satellite), and Ld^ij is a delay caused 
by clock estimate errors at both the satellite and the re- 
ceiver. All other symbols are the same as in Section II, 
and the ionospheric delays have been assumed to be cali- 
brated. Similar to Eqs. (la) and (lb), Eqs. (8a) and (8b) 
can be parameterized by setting {Ltr^s) = Ag Ltr,\,z 3-nd 
— Ai Ltrj^zi where ^re the statistically av- 

eraged tropospheric delays in the zenith direction at jth 
site. By assuming that the observable errors have zero 
means and a variance-covariance matrix IV” \ the param- 
eterized Eqs. (8a) and (8b) can be solved by least-squares 
estimation for the geometrical delays Lg g and Lg^ij (or 
rather, by using geometrical delay modeling for parameters 
which determine Lg^g and j j), clock errors Lcj^ijy and 
zenith delays Ltrj^z^ As was discussed in Section II, it is 
this parameterization, accompanied by the statistical eval- 
uation of W~^y that allows the solution to produce best 
estimates. Note also that once Lg^g has been estimated, 
the best estimate for line-of-sight can in principle be 
obtained as the difference Ltr,s = ^tot,s ~ 

For real data, the accuracy of path delay estimates, in 
addition to tropospheric inhomogeneities, will be affected 
by other error sources, including instrument and multi- 
pathing noise and uncertainties in geometrical delay mod- 
eling. Neglecting the specifics of geometrical delay model- 
ing and error statistics, the effect of an overall uncertainty 
level contributed by various error sources has been mod- 
eled by adding to each cov (Ltr^iy Ltr,i) a term assumed to 
be uncorrelated between different directions and propor- 
tional to the air mass. Figure 8 shows the effect of the 
assumed 0.1-, 0.3-, and 1-cm uncertainty levels.^® Note 
that these levels exceed the tropospheric inhomogeneity- 
induced errors at Eg = 50, 22, and 12 deg, respectively. 
That is, to achieve the tropospheric inhomogeneity-limited 
accuracy at some elevation, the additional uncertainty 
must be reduced below the inhomogeneity-induced error 
at that elevation. The dashed curves in Fig. 8 show the 
error resulting from the use of Cj’s which were optimized for 
zero measurement and modeling uncertainty. The errors 
are bigger than for the optimized c,’s by about 15 percent, 


Analysis of recent GPS data shows that assumption of uncorre- 
lated errors is not entirely correct. Efforts to understand the error 
statistics are xmder way. 


V. Conclusions and Recommendations 

This article has outlined a mathematical procedure 
which has the potential for minimizing the tropospheric 
inhomogeneity-induced mapping error for estimated path 
delays along arbitrary lines of sight by using observables 
from different directions. The main principle of the proce- 
dure is the application of least-squares estimation to the 
combined set of observables in the calibration directions 
and the direction to be calibrated, and the use of an ob- 
servable variance- covariance matrix during estimation. A 
numerical example was given for a set of calibration direc- 
tions assumed to coincide with directions from a ground- 
based receiver to GPS satellites visible at Goldstone. For 
these directions, assuming an error-free calibration instru- 
ment and zero geometrical delay mismodeling, and using 
the Kolmogorov turbulence model, the azimuthally aver- 
aged error is found to be about 1 mm in the elevation 
range from about 40 deg to zenith. At elevations less than 
40 deg, the error increases with decreasing elevation, reach- 
ing about 1.2 cm at 10 deg. Because of its stochastic ori- 
gin, the inhomogeneity-induced error cannot be removed 
by improved modeling or instrument design, and repre- 
sents the ultimate accuracy for line-of-sight estimates. 

The minimum error was compared to the error obtained 
by using zenith delays which were averaged over many cal- 
ibration directions, and to error obtained by using line-of- 
sight estimates optimized for zero wind. The zenith map- 
ping error is nearly twice as large (in the absence of other 
error sources) as the minimum error, indicating that the 
effort spent to minimize the error is worthwhile. On the 
other hand, the error increase resulting from the use of zero 
wind was less than 30 percent, indicating that it may not 
be necessary to include the wind in the estimation proce- 
dure. Signal integration reduces the error; the integration 
time required to halve the instantaneous error is less than 
8 min for elevations between 30 deg and zenith, reaching 
about 12 min at 10 deg. The use of multiple on-site re- 
ceivers will not help to reduce the error at low elevations 
where it would be needed the most. 

The suggested procedure could be used to estimate path 
delays in the direction of the DSN antenna by using GPS 
data. The present strategy for GPS data analyses en- 
tails construction of a tropospheric zenith delay by av- 
eraging over all satellites; the estimated zenith delay error 
is at the 1-cm level due to the combined effect of tropo- 
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spheric inhomogeneities, uncertainties in geometrical de- 
lay modeling, and measurement (thermal noise and mul- 
tipathing) uncertainties [9]. The combined effect of ad- 
ditional uncertainties has been modeled by adding to the 
tropospheric variances a term assumed to be uncorrelated 
between various measurements. As expected, the accuracy 
of the estimate was limited by this additional term when- 
ever it exceeded the tropospheric inhomogeneity-induced 
level. Further work should be performed to establish levels 
and statistics of additional error sources, incorporate their 


correlations into the analysis, and, if possible, reduce their 
effect on path delay estimates. 

In spacecraft tracking, observing epochs and scan 
lengths are specified by mission considerations, and the 
tracked directions differ from satellite directions. How- 
ever, experiments could be designed in which the direction 
to be calibrated coincides with one of the satellite direc- 
tions. For example, by directing WVRs towards satellites, 
the WVRs themselves could be calibrated. 
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Fig. 1. The coordinate system for DSN and calibration directions. 



Fig. 2. The instantaneous, tropospheric fluctuation-induced path delay error versus the observed elevation, 
for (a) /? = 0 m and (b) R = 1000 m. The error range Is due to the 360-deg range in the observed azimuth, 
the azimuth of R, and satellite constellations. The dashed curves are the azimuthaliy averaged error. 




PATH DELAY ERROR, cm 



Fig. 3. The tropospheric fluctuation*!nduced error for one cal- 
ibration line pointed In the (f>f = 45-deg, £/ = 28.4-deg direction 
(A? = 0, error-free calibration Instrument). The error with optimized 
c/’s (solid and dashed-line curves) Is less than the zenith mapping 
error (dash-dot curve) by an amount that depends on the eleva- 
tion. 



Fig. 4. Tropospheric fluctuation-induced error and the error range 
(averaged over four principal wind directions) versus the signal 
Integration time T. The azimuth of ^ (A? = 200 m) Is along the 
azimuth of <f>s, v = 10 m/sec, and the satellite constellation is for 
12:00 a.m., July 1, 1992; C is the turbulence strength constant. 
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Fig. 5. Tropospheric fluctuation-induced error with c/’s optimized 
for three calibration instruments positioned at corners of an equi- 
lateral triangle centered at f? = 200 m around the DSN antenna 
(dash-dot line); one instrument located at R = 0 (solid line); and 
one instrument at 200 m (dashed line). 



Fig. 6. Tropospheric fluctuation-induced, azimuthally averaged 
rate error versus scan duration {R = 200 m, v = 10 m/sec). The 
error caicuiated by using the optimized c/'s becomes smailer than 
the uncalibrated error and the error for zenith mapping at a Tsc 
that depends on elevation. 
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DELAY RATE ERROR, sec/sec 



ELEVATION, deg 

Fig. 7. Tropospheric fluctuation-induced rate error for scan dura- 
tions Ts - 100 sec and 1000 sec versus £5 when = 0. The satel- 
lite constellation used (12:00 a.m., July 1, 1992) has one satellite 
in the <j>i - 45-deg direction. 



ELEVATION, deg 


Fig. 8 . The combined effect of additional uhcdrrelifted uncertainty^ 
levels on the total (azimuthally averaged) error f or R = 200 m. The 
solid-line curves were computed by using c/s optimized for the 
additional uncertainty; the dashed lines used c/s optimized for 
zero additional uncertainty. 
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Appendix A 

Tropospheric Correiations Using the Koimogorov Turbuience Model 


Path delay correlations are evaluated by writing the path delay at site i and time t as the integral 

Ltr,i{t) = / X{ri,t)dri (A-1) 

Jo 

where the air mass Ai = l/sinE'j, hy is the wet troposphere slab height, and index of refraction — 1 at 

location and time t. For example, the correlation 


cov -\-T)) 


+ "F)) — {Ltr^i) {Lfrj)) 


(A-2) 


is evaluated by substituting Eq. (A-l) and the expression (Eq. (A. 3) of [1]): 

{x{ri,i)x{rj,t + T)) = (x^) - - fj + v T) 


(A-3) 


where vis the wind velocity, and D^{fi -fj+vT) = ^(x(ri,<) - + 7’))^) ‘s the structure function for inhomo- 

geneities correlated both spatially and temporally. By interchanging the order of integration and ensemble averaging 
and setting dri ~ A{dz and dvj = Ajdz\ Eq. (A-2) becomes 


cov {Ltr,t{t),Ltrj{i -hT)) 



D^{fi - Tj -hu T) 


(A-4) 


where the variance = (x^) - (x)“ is independent of spatial coordinates and is obtained by letting the distance R go 
to infinity in 


Assuming that the troposphere is described by the Kolmogorov turbulence model, the structure function for the frozen 
inhomogeneities is 


D^(^R -|- = 


[x{r,t) - x{ 


r “f R-\- V t 


»’) 


l + {\R-^i7t\/L,)V^ 


(A-6) 


where R is the spatial interval over which the structure function is evaluated, the saturation scale length Lg 3000 
km, and the turbulence strength C = 1.1 x [1,8] (for the wet'slab height hy czi 2 km, corresponding to about 

a 6-cm zenith wet path delay at Goldstone). 

Equations of the same type were also used to quantify the effect of dry fluctuations. Assuming that the dry 
fluctuation is one-third of the wet fluctuation [1] with the scale height — 8 km, the dry turbulence strength 
Cd {Cy /3){hy / — 9.2 x 10”^. For uncorrelated fluctuations, dry and wet errors add in quadrature. 
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Appendix B 
The Delay Rate Error 


The delay rate error was evaluated as the square root 
of the variance of the estimated delay rate, ^(T^c)? which 
is the slope of the path delay over the scan duration T,c 
[1]. By dividing T,c into N sections of equal length, and 
assuming that the time origin is in the middle of T,e, a 
linear fit to points at the centers of the sections gives the 
following formula: 


where Ltr,s{^), and Ltr,s are the actual 

and estimated tropospheric path delays in the observed, 
{Esy(f)s)^ direction, at times t and o, respectively. 

By using Ltr^s = J2k ^kLtr,k (see the discussion follow- 
ing Eq. (5) of the main text), Eq. (B-2) was evaluated by 
substituting Eq. (B-3) into Eq. (B-2). This yields 


5(T,c) = (B-1) 

where L{ii) is the delay at the time point tj at the center 
of the zth section. By using a similar expression for the 
actual delay rate, B{Tsc)y tbe delay rate variance is derived 
as 


^\{Tsc) — (^.^?)2 {Ltr,s{U 


+ ^ ^ CkCl COV {Ltr^kiU - ij), Ltrj) 
h I 

(cOV g (/ j ) ) Eif^k) 

k 


1 




(B-2) 


4- COV (Ltr,k(ti - tj 



(B-4) 


where i, j = 1, N, and the path delay variance, tr? 
is 


where Ttr,Jfe(0 tropospheric path delays in (Ek,(f>k) di- 
rections at time i, Eq. (B-4) was evaluated for i =!,. •> 
3 equally spaced sections within Tgc (ti = —Tgc/^i to = 
0, t 2 = T,c/3) by using the Kolmogorov turbulence model 
described in Appendix A. 
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